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Anisotropic magnetoresistance (AMR) is a relativistic magnetotransport phenomenon arising from 
combined effects of spin-orbit coupling and broken symmetry of a ferromagnetically ordered state 
of the system. In this work we focus on one realization of the AMR in which spin-orbit coupling 
enters via specific spin-textures on the carrier Fermi surfaces and ferromagnetism via elastic scat- 
tering of carriers from polarized magnetic impurities. We report detailed heuristic examination, 
using model spin-orbit coupled systems, of the emergence of positive AMR (maximum resistivity for 
magnetization along current), negative AMR (minimum resistivity for magnetization along current), 
and of the crystalline AMR (resistivity depends on the absolute orientation of the magnetization 
and current vectors with respect to the crystal axes) components. We emphasize potential qualita- 
tive differences between pure magnetic and combined electro-magnetic impurity potentials, between 
short-range and long-range impurities, and between spin-1/2 and higher spin-state carriers. Con- 
clusions based on our heuristic analysis are supported by exact solutions to the integral form of the 
Boltzmann transport equation in archetypical two-dimensional electron systems with Rashba and 
Dresselhaus spin-orbit interactions and in the three-dimensional spherical Kohn-Littinger model. 
We include comments on the relation of our microscopic calculations to standard phenomenology 
of the full angular dependence of the AMR, and on the relevance of our study to realistic, two- 
dimensional conduction-band carrier systems and to anisotropic transport in the valence band of 
diluted magnetic semiconductors. 



PACS numbers: 72.10.-d, 72.20.My 

I. INTRODUCTION 

Advanced theoretical approaches and experiments in 
new unconventional ferromagnets have recently led to a 
renewed interest in the relativistic, extraordinary magne- 
totransport effects. There are two distinct extraordinary 
magnetoresistance coefficients, the anomalous Hall effect 
(AHE) and the anisotropic magnetoresistance (AMR). 
The AHE is the antisymmetric transverse magnetoresis- 
tance coefficient obeying pxyCM.) = —pxy{—'M.), where 
the magnetization vector M is pointing perpendicular 
to the X, y plane of a Hall bar sample. The AMR 
is the symmetric coefficient with the longitudinal and 
transverse resistivities obeying, p^xi]^) = M) and 

Pa;j,(M) = pxy{—'M), where M has an arbitrary orien- 
tation but in most studies it lies in the x — y plane. 
Numerous works have explored the origins of the AHE; 
for reviews see e.g. Refs. Diluted magnetic 

semiconductors became one of the favorable test bed 
systems for AHE investigation^'^'^i^i^ii^ due to their 
tunability and the relatively simple, yet strongly spin- 
orbit coupled Fermi surfaces i ^^'^^ An even more system- 
atic and comprehensive understanding of the AHE on 
a model level has been obtained by considering two- 
dimensional semiconductor systems with archetypical 
spin-orbit interactions (SOIs) of the Rashba and Dres- 



selhaus type4ii^^ii^^iia8i9^^^^ 

Despite the long history and importance in magnetic 
recording technologies, the AMR has been studied less 
extcnsively,^i^i2^i^i2^ Similar to the AHE, it has been 
recently argued that the analysis of the AMR can be 
significantly simplified in diluted magnetic semiconduc- 
tors hkc (Ga,Mn)As.^'^ Two distinct microscopic mech- 
anisms have been identified that can lead to anisotropic 
carrier life-times in these systems: One combines the 
spin-orbit coupling in an unpolarized carrier band with 
scattering off polarized magnetic impurities while the 
other emphasizes polarization of the carrier band itself 
and does not require magnetic nature of the scatterers. 
(Note that apart from life-times, the AMR may also arise 
from anisotropic group velocities4^) Although acting si- 
multaneously in real systems, theoretically both mecha- 
nisms can be turned on and off independently and it was 
found3^ that the scattering of spin-orbit coupled band 
carriers from magnetically polarized impurities should 
dominate in the diluted magnetic semiconductors. Build- 
ing on the analogy with AHE studies we seek further 
insight into the basic physics of this AMR mechanism 
by focusing on the archetypical spin-orbit coupled two- 
dimensional systems. 

Using the relaxation-time approximation (RTA) and 
starting with the Rashba and Dresselhaus models we 
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show in Sec. [TTlhow the sign of the AMR can be inferred 
by inspecting the spin texture of the spin-orbit coupled 
Fermi surface. We point out that impurities containing 
polarized magnetic potential only or containing a com- 
bined electro-magnetic potential can yield distinct AMR 
phenomenologies. Examination of the Rashba and Drcs- 
selhaus models allows us to draw separate links between 
the spin-texture and the non-crystalline and crystalline 
AMR components where the non-crystalline AMR de- 
pends on the relative angle between M and current I 
while the crystalline AMR has an additional dependence 
on the absolute orientation of M and I in the coordinate 
system of the crystal axes. We conclude the qualita- 
tive discussion in Sec. Hi] by illustrating in the Rashba- 
Dressclhaus system a potentially important effect on 
AMR of long-range impurities, and in a spherical Kohn- 
Luttinger model^- the effect of carriers with higher spin 
state. Analysis of these effects relates our work to pre- 
vious theoretical studies of the AMR in (Ga,Mn)As di- 
luted magnetic semiconductors j^i^i^i^^ The validity of 
the heuristic analysis of the AMR is confirmed in Sec. IIIII 
where we explain the relation between the RTA and 
the exact solution to the integral Boltzmann cquationf^ 
Quantitative results for the AMR arc derived in this Sec- 
tion and Appendix for the Rashba model and for the 
Dresselhaus model with short-range electro-magnetic im- 
purities and for the combined Rashba-Dresselhaus model 
with arbitrary strength of the two SOI terms and with 
short-range magnetic impurities. In Sec. lIVI wc comment 
on the relevance of our model calculations to realistic 
two-dimensional semiconductor structures. 



II. HEURISTIC LINK BETWEEN SPIN 
TEXTURES AND IMPURITY POTENTIALS 
AND THE AMR 

We limit our discussion in this section to AMRs defined 
as the relative difference between longitudinal resistivi- 
ties for magnetization aligned parallel and perpendicular 
to the current direction. In situations discussed below, 
the transverse resistivity vanishes and we can define 



AMR: 



ipl+pf)/2 ial+a^)/2 



(1) 



where (cj) and p-^iu-^) is the longitudinal resistiv- 
ity (conductivity) for M || I and for M _L I, respec- 
tively, and the subscript I labels the orientation of cur- 
rent with respect to crystal axes. (The relation of our 
microscopic theory to the standard phenomenology of 
the full angular dependence of the AMR will be com- 
mented upon in Sec. IIIII ) Our heuristic analysis of the 
AMR defined in Eq. ([T|) is based on the RTA and on as- 
suming a proportionality between resistivity and the 1st 
order Born approximation to elastic scattering probabil- 
ities from the state with the group velocity along I. Fur- 
thermore we consider only the strongest contribution to 



the transport life-time which comes from back-scattering, 
i.e., from transitions into states with group velocity op- 
posite to I. We use these approximations and consider 
several archetypical spin-orbit coupled Fermi surfaces to 
elucidate the relation of the spin-texture and nature of 
the impurity potential to various fundamental aspects of 
the AMR phenomenology. 



A. AMR in the Rashba model 

We start with the two-dimensional electron system 
with Rashba SOI which yields positive AMR indepen- 
dent of the current orientation in the crystal, and demon- 
strate the potential qualitative difference between pure 
magnetic short-range impurity potential oc gm ■ s/s and 
a combined electro- magnetic potential cx 1 + e^j - s/s. 
Here denotes the magnetization unit vector and s is 
the carrier spin operator. For electrons with s = 1/2, 
the operator s/s can be represented by the 2x2 Pauli 
matrices cr — {a^, (Jy,(Jz)- 

The tangential spin-texture along the Fermi contour of 
the Rashba Hamiltonian, 



2l2 



Ht 



2m 



(2) 



is shown in Fig. [ija). The spinors on the majority (— ) 
and minority (-I-) Rashba band are given by |k±) = 
(1, =F*e*^); where ta,n0 = ky/kx- From now on the coordi- 
nate system is chosen in such a way that x, y, and z direc- 
tions coincide with [100], [010], and [001] crystallographic 
axes respectively, as shown in Fig. [2l Assuming current 
along i-dircction, we can infer the back-scattering am- 
plitudes of the states with the group-velocity (k-vector) 
parallel to the current by recalling the following proper- 
ties of the scattering matrix elements: 



U l^^x|i) = o 

a Wv\ i) = 1 



(T \<y.\ i) = 1 
(T 1^,1 i>=o. 



(3) 



Here we labeled the spinors by arrows whose orienta- 
tion can be directly compared to the spin-textures de- 
picted in Fig. [Ija). The allowed back-scattering pro- 
cesses, according to the relations in ([3]), are highlighted 
in Fig. HKb) for the pure magnetic impurity potential. 
When magnetization points along the s-direction (i. e. 
to the right in Fig. [Hb)), Bm ■ cr = cTx and the back- 
scattering of states moving along the i-direction is due 
to majority-to-majority and minority-to-minority band 
transitions. In the case of magnetization parallel to 
the y-direction, eu ' ^ = '^y and back-scattering is due 
to majority-to-minority and minority-to-majority transi- 
tions. In the hmit of k- » fc+, these figures suggest 
that back-scattering is strongly suppressed for M _L I 
implying low resistivity in this configuration compared 
to the M II I case. The AMR defined in Eq. ^ is there- 
fore expected to have positive sign in the Rashba model. 
Quantitative Boltzmann equation calculations presented 



3 




FIG. 1: Rashba model and (a) its spin texture along the Fermi 
contours. Dominant scattering channels for the states with 
group velocity pointing to the right when (b) magnetic and (c) 
electro- magnetic impurities (see text) constitute the prevalent 
source of momentum relaxation. Note the indicated directions 
of impurity polarization. The current flow is directed from left 
to right. The reader might consider the limit k- 3> k+ for 
better understanding of this and subsequent figures. 



in Sec. IIIII confirm the positive AMR for all fc_ > fc+. 
They also confirm the vanishing magnitude of the AMR 
in the weak SOI, large Fermi energy limit w fc_) 
which is discerned directly from our pictorial represen- 
tation of the allowed backs-scattering transitions consid- 
ering nearly degenerate majority and minority Rashba 
bands in Fig. [T]^b). 

The behavior of AMR in the limit of degenerate 
Rashba bands, while keeping the tangential spin tex- 
tures, is qualitatively altered when the impurity potential 
contains magnetic and non-magnetic components (e.g. 
for Mn acceptors in III-V semiconductors). Replacing 
ax,y with 1 + crx,y in the relations ([3]) allows us to illus- 
trate this by again considering the transitions that con- 
tribute to the back-scattering; note that this does not 
describe the situation where there are two distinct types 



of impurities^ (such as phonons and charge-neutral mag- 
netic impurities). As highlighted in Fig.[T]Jc), there is now 
always one of the Rashba bands in which back-scattering 
is absent for M _L I, independent of the difference be- 
tween and k^. For M || I, back-scattering occurs 
in both bands and each of the states moving along the 
current can scatter to both majority and minority band 
states. This implies large positive AMR even in the limit 
of fc+ w fc_ . 

Finally we point out that the circular symmetry of the 
Rashba spin-texture makes the model a prototype real- 
ization of a purely non-crystalline AMR system. The 
AMR is independent of the orientation of current in the 
coordinate system of crystallographic axes and depends 
only on the relative angle between M and I. 




FIG. 2: Dresselhaus model and (a) its spin texture. In order 
to determine the current and the AMR along the [100] and 
[110] crystallographic directions we focus on the states with 
group velocities pointing in the respective directions, (b) and 
(c). Dominant momentum relaxation channels for these states 
and scattering on magnetic impurities are indicated on the 
bottom panels. 



B. AMR in the Dresselhaus model 

The tangential spin- 1/2 texture of the Rashba model 
represents arguably the simplest host for a positive 
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purely non-crystalline AMR. The Dresselhaus SOI can 
be viewed as a minimal model demonstrating the link 
between a radial spin-1/2 texture and a negative AMR, 
and illustrating the emergence of crystalline AMR. The 
Dresselhaus Hamiltonian, 

Hd = + Picr^kj; - Oyky) , (4) 

yields the majority and minority cigcnstates, |k-j-) = 
(l,±e~'^), whose spin orientations along the respective 
Fermi contours are depicted in Fig. HJa). We can use 
the same analysis of the back-scattering amplitudes as in 
the previous subsection to link this spin texture to the 
expected basic AMR phenomenology in the Dresselhaus 
model. 

In Fig. mb), we consider the case of current flowing 
along the £-direction ([100] crystal axis) and scattering 
from impurities carrying the short-range magnetic poten- 
tial only. Using the same representation of the spinors as 
in Eqs. ([3]) we can write 

(^|a,|^) = l {^K\^)=0 

This implies that for magnetization parallel to the cur- 
rent direction, back-scattering is due to majority-to- 
minority and minority-to-majority band transitions while 
for magnetization perpendicular to the current, allowed 
transitions are the majority-to-majority and minority-to- 
minority. The low-resistivity and high-resistivity magne- 
tization orientations therefore switched places compared 
to the Rashba model and the AMR becomes negative. 

The spin-texture of the Dresselhaus model is not cir- 
cularly symmetric, however. It evolves from radial for 
k parallel to the [100] or [010] crystal axes to tangential 
for k parallel to the [110] or [110] diagonals, as shown 
in Fig. UK a). The back-scattering amplitudes for cur- 
rent along the diagonal, highlighted in Fig. ^c), are 
hence identical as in the Rashba model, implying pos- 
itive AMR for this current direction. The lower sym- 
metry of the Dresselhaus model does not give rise to 
anisotropy in the resistivity of the system in the ab- 
sence of magnetization^^ However, when magnetiza- 
tion is present the system acquires a crystalline AMR 
which reflects the underlying cubic symmetry of the spin- 
texture. We remark that both the negative and posi- 
tive AMRs of the Dresselhaus model vanish in the limit 
of fc+ w Also in analogy with the behavior of the 
Rashba model, the AMRs with the respective signs are 
recovered in this limit when the pure magnetic impurity 
potential is replaced with the combined electro-magnetic 
potential (see Sec. IIIII and Fig. (5]). 

C. AMR in the Rashba-Dresselhaus model with 

I"! = 1/31 

We now briefly comment on the potential importance 
of long-range nature of the impurity potential on the ba- 



(a) 



s=l/2 




FIG. 3: (a) Spin texture along the Fermi contours of Rashba- 
Dresselhaus model with a = /3. The AMR is zero for any type 
of short-range impurities. However, for long-range magnetic 
impurities the scattering amplitudes depend on the momen- 
tum transfer (illustrated by the length of the arrows) and 
non-zero AMR arises for current both along (b) [110] and (c) 
[110] crystallographic directions. 



sic AMR phenomenology. For the demonstration of this 
effect, a singular model combining Rashba and Dres- 
selhaus SOIs with |q;| = \(3\ is particularly suitable. 
The Hamiltonian containing Rashba and Dresselhaus 
spin-orbit coupling terms of equal strength has singular 
properties^^i^ii^ (in particular additional symmetries). 
The internal spin-orbit coupling field has a k-vector in- 
dependent orientation (along the [110]-axis for « = /?). 
Spins on one circular Fermi contour are aligned parallel 
to this field while on the other contour they take the anti- 
parallel alignment. Additionally, as shown in Fig. [3Ja), 
this singular SOI shifts the two equal-size Fermi contours 
with respect to each other along a direction perpendicu- 
lar to the direction of the spin-orbit field. 

Because of the rigid spin-texture of the |a| = |/?| 
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Rashba-Drcsselhaus model on two mutually shifted but 
otherwise identical circular Fermi contours, the back- 
scattering amplitudes for a short-range impurity poten- 
tial are independent of both the relative angle between 
M and the group velocity of the state moving along I, 
and of the direction of current with respect to crystal 
axes. The AMR therefore completely vanishes in this 
model. Nevertheless, Figs. |31Ib),(c) illustrate that the 
AMR, including its crystalline component, is recovered 
when the scattering amplitudes pick up a dependence on 
the transferred momentum, i.e., for impurities carrying a 
long-range electro-magnetic potential. 



D. AMR in the spherical Kohn-Luttinger model 

We conclude our excursion into the basic phenomenol- 
ogy of AMR, produced by scattering of spin-orbit coupled 
carriers from polarized magnetic impurities, by consider- 
ing higher spin state of the carriers. We show that seem- 
ingly identical spin-textures can result in opposite sign of 
the AMR for spin-1/2 and higher spin carriers, and argue 
that the AMR can have opposite sign when carriers with 
higher spin are scattered from a pure magnetic or from a 
combined electro-magnetic potential. Again seeking the 
minimal SOI model on which this AMR phenomenology 
can be demonstrated without performing detailed trans- 
port calculations we choose the four-band spherical three- 
dimensional Kohn-Luttinger Hamiltonian for total angu- 
lar momentum j = 3/2 carriers. 



KL 



2m 



(7i + T^72)fc'- 272(k-j)2 



+ hj, (6) 



with h 0. The k^, ky plane (with infinitesimal fc^) spin- 
textures depicted in Fig. [DJa) are obtained by realizing 
that the spin operator s = j/3 in the four-band model, by 
defining the momentum quantization axes parallel to k, 
and considering only the jk = ±3/2 bands (heavy holes). 
The infinitesimal exchange field h in Eq. ^ is included 
to lift the degeneracy of these two bands, and 71 and 72 
are the Luttinger parameters specific to the particular 
semiconductor valence bands for which Hkl is derived 
from the conventional k • p approximationi22i2& 

Unlike the spin-1/2 Dresselhaus model, the radial spin 
texture in the j = 3/2 Kohn-Luttinger model yields a 
positive AMR for purely magnetic scatterers. This can 
be illustrated using an analogous representation as in 
Eqs. ([5]) to relate the scattering amplitudes for impurity 
potential oc gm -8/8 = gm • j/j and the spin-texture. For 
the j = 3/2 carriers we obtain^! 



^0 
= 



\jy\ 



= 
= 



(7) 



This implies, as highlighted in Fig.UJb), that for magne- 
tization parallel to the current direction, back-scattering 
is due to majority-to-minority and minority-to-majority 




J = 3/2 




FIG. 4: (a) Cross-section (parallel to the 
the 3D radial spin texture belonging to the 
bands of the Kohn-Luttinger Hamiltonian. 
surfaces are sketched with different sizes for 
the Hamiltonian ((Gjl implies k- — k+ as h ^ 
scattering channels for magnetic impurities 
ence to Fig. [IJb). (c) The same as (b) for 
scatterers. 



kx,ky plane) of 
two lower-energy 
The two Fermi 
clarity, although 
0. (b) Dominant 
note the differ- 
electro-magnetic 



band transitions as in the case of spin-1/2 carriers. How- 
ever, for magnetization perpendicular to the current, 
there are no allowed back-scattering transitions in con- 
trast to the spin-1/2 Dresselhaus model in Fig. [2Kb). This 
makes now the latter configuration the low-resistivity 
state and AMR for the radial spin-texture of the Kohn- 
Luttinger model becomes positive for pure magnetic im- 
purity potential even for k+ ~ k-. Boltzmann equation 
calculation of the AMR presented in Appendix [B] (and 
also an independent calculation based on the Green's 
function formalism^^) again confirm our heuristic con- 
clusion of Fig. IHb). 

On the other hand, electro-magnetic scatterers oc 
|l + ix,y produce negative AMR in the Kohn-Luttinger 
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niodel2ii2^ in the very same way as it is shown in Fig. [6Kb) 
for the Dresselhaus model, and in both cases, this be- 
havior can again be inferred using relations ([7|) and ([5|) 
with jx^y and ax,y replaced by |l + jx,y and 1 + (7x,y, 
respectively. Dominant scattering channels which sug- 
gest that AMR < are summarized in Fig. HJc). Con- 
trary to the Dresselhaus model ^ , the SOI of the Kohn- 
Luttinger model ^ in combination with polarized scat- 
terers therefore can produce AMR of either sign, e.g. 
depending on the carrier-density-controlled screening of 
the impurities i^i^ This qualitative difference between 
Dresselhaus and Kohn-Luttinger models highlights the 
fact that knowledge of spin textures, such as Figs, ^a.) 
orSJa), may not be sufficient to analyze the scattering 
properties of the model and appropriate matrix elements 
such as Eqs. (O or ([7]) should always be verified. 

III. QUANTITATIVE RESULTS FOR THE AMR 
IN THE RASHBA-DRESSELHAUS MODEL 

The AMR analysis in the previous Section utilizes the 
RTA (in fact only the back-scattering term of the RTA) 
which, in general, is not a rigorous theory approach for 



where v^^k = dEi^k/dhk is the group velocity, /(i,k) is 
the non-equilibrium distribution function, and i = ± is 
the band index. The transition probabilities in the 1st 
order Born approximation arc given by 

«;(^,k;^^k') = ^|(^,k|y|^',k')|^ (9) 

where V is the strength of the short-range scattering po- 
tential of impurities with density n. Energy conservation 
during elastic scattering processes was already incorpo- 
rated into the right hand-side of Eq. ([S]). 

In the Rashba model, J^i' ''^ih a constan-fc^ 

for a short-range electric potential, oc 1, or magnetic 
potential, V oc cm ■ cr- In the limit of nearly degenerate 
bands, Ei ]^ fa Eii we can find a solution of Eq. ([S]) in 
the RTA form, 

/(z, k) - /0(£;,k) = c|e|E . v.k^^^J^ . (10) 



anisotropic systemsi^ It is therefore desirable to calcu- 
late the AMR beyond the RTA, not only to obtain quan- 
titative predictions but also to confirm the validity of the 
basic AMR phenomenology inferred above. As in Sec. [Til 
we will employ the 1st order Born approximation for cal- 
culating the scattering probabilities but will solve the 
corresponding integral Boltzmann equation exactly. To 
provide better physical insight we start with explaining 
the relation between the RTA and the full semiclassical 
Boltzmann theory for the two-dimensional SOI systems. 
Exact analytical solutions to the Boltzmann equation 
are then derived for Rashba and Dresselhaus model with 
short range electro-magnetic impurity potentials and for 
the combined Rashba-Dresselhaus model with arbitrary 
a and (3 and with magnetic impurities. 



A. Relation between RTA and integral Boltzmann 
equation in the Rashba model 

Because the equilibrium Fermi distribution f^{Eix) is 
a function only of energy, we can write the Boltzmann 
equation^ in d = 2 dimensions as 



(8) 

I 

Plugged in Eq. ([5]), the second term on the right- 
hand side drops out because of the independence of 
w(i, k; i', k') on k' and because the group velocity 
averages to zero over the Fermi contour, and the first 
term gives 

\^ j 7S^E^(*'i^;*''W£;,'.k'-i?.,k)^i. (11) 

i' 

The electrical current within the semiclassical linear re- 
sponse, given by 

j = E / ^^^'^ - 1^)] ' (12) 

is exactly proportional to the quasiparticle broadening 
life-time t in this case. Same RTA form of the Boltz- 
mann equation applies also to the Rashba-Dresselhaus 
model with |q;| = \(3\ because the rigid spin-texture of 
this singular case implies constant transition probabili- 
ties for any short-range electro-magnetic potential. 



-|e|E.v,,k^^^^^ = - I -0j^u;(*,k;i', k')<5(i?,,,k'- i?«,k)[/(*,k) -/(*', k')] 
= ~[/(z,k)-/°(£;,,k)]5^ j -^wi^,k■,^'X)SiE,,^,-E,^^) 
+ J ^T.'^{hk;t',k')6{Ee^u'~E.,,^)[f{t',-k')~f{E,,,^,)] 
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In the Rashba model with non-degenerate bands, 
Ei.k ^ Ei^^]^, the RTA solution pI7|) to the Boltzmann 
equation can still be found for a non-magnetic potential, 
V (xt. The scattering probability w(i,k;i',k') depends 
in this case on the magnitude of the transition angle, 
\9 — 6'\. It implies that from the product, 

E-v,/ k' = ViE-Vk— cos(9-0')+v,(zxE,)-Vk— sui(9-9') 

(13) 

the transverse term oc sm{9 — 9') does not contribute to 
the second term on the right-hand side of Eq. The 
longitudinal term oc cos{9 — 9') contributes to Eq. ^ and 
the Boltzmann equation takes a modified RTA form with 

I ^ /^E^(*'k;z',k')<5(i?,',k'-ii;.,k) 



bmce also J^"" d9'wii,9;i',9') is a constant independent 
of 9, the first term on the right-hand side of Eq. ^ im- 
plies that /(j, k) — f^{Ei,\f^) must contain term E • Vi,fc(6') 
and the second term on the right-hand side of Eq. ([5]) 
implies that /(i, k) — /"^(i^i.k) must contain harmonics of 
w{i, 9; i', 9') which in both cases happen to be just cos 61 
and sin 9. No higher order Fourier components can con- 
tribute to the non-equilibrium distribution function in 
this case and Eq. ([S]) can be solved analytically. 

The AMR of the Rashba model with magnetic impu- 
rity potential is summarized in the first column of Tab. |T] 
and also plotted in Fig. [5] as a function of the ratio 
Epfi^ /{ma^). Here Ep = corresponds to the minority 
Rashba band being just depleted and Epti? /{ma^) » 1 
to nearly degenerate i = ± Rashba bands. Consistent 
with the qualitative results of Sec. [ll] we find a positive 
AMR which vanishes as the radii of the minority and 
majority band Fermi contours approach each other. 



1 



■ cos W — 



1 

Ttr 



(14) 



Electrical current is now proportional to the transport 
life-time which gives larger weight to larger angle scat- 
tering transitions. 

The transport life-time form of the Boltzmann equa- 
tion has been the basis of qualitative discussions in Sec.HIl 
where wc further simplified the analysis by considering 
only the leading contribution to current in Eq. (|12p from 
states with v.^^k || E. For all spin-textures and orienta- 
tions of E and M considered in Sec. |TT1 w{i, k; z', k') de- 
pends only on |0 — 0'| for the special h-states with group 
velocity parallel to the electric field. This justifies the 
internal consistency of the RTA based analyses in Sec. HIl 
and explains their qualitative validity. 



B. Solution to the Boltzmann equation for the 
Rashba-Dresselhaus model 

To obtain quantitative AMR predictions we need to 
perform the full k-space integration in the expression 
(fT2|) for the electrical current. For arbitrary k- vector 
and other than the few special cases discussed in the pre- 
vious subsection (which all happen to give zero AMR), 
the integral of the transverse term in Eq. may not 
vanish and/or the integrated scattering probability in the 
first term on the right-hand side of Eq. ([8]) may not be 
independent of k. In these cases the RTA form of the so- 
lution to the Boltzmann equation fails. For the Rashba- 
Dresselhaus model we can, nevertheless, find the exact 
solution to the Boltzmann equation in an analytic form 
which allows us to directly compare the corresponding 
quantitative AMR predictions with the qualitative re- 
sults of Sec. n 

The method has been previously derived^ for pure 
Rashba model in which the angular dependence of the 
scattering probability function for the short range mag- 
netic potential, e.g. V oc ax, is given by 
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FIG. 5: Pure Rashba system with magnetic impurity, AMR 
as a function of tlie Fermi energy Ep in units of ma^ /H^. 

For Rashba model with the electro-magnetic poten- 
tial, e.g. V cc t + ax, the integral J^^ d9'w{i, 9; i' , 9') oc 
1 -|- ii' sin 9 is not a constant which implies the presence 
of higher order Fourier components in /(z, k) — f^{Ei ]f^). 
Still an analytical form can be found for the distribution 
function, see the note added in proof of Ref. [H- Analo- 
gous arguments apply also to the Dresselhaus model with 
electro-magnetic impurities. The dependence of AMRs 
in the two models as a function of the ratio a of the 
electrical and magnetic parts of the impurity potential 
V al + Cm • cr in the limit of nearly degenerate bands 
and for current along the [100]-axis is given by 



AMR 



±2a^ for \a\ < 1 
±2/a2, for \a\ > 1 



(16) 



■w{i, 9; i' ,9') (X 1 ~ ii' {cos t 



(15) 



where +/— corresponds to the Rashba/Dresselhaus 
model. For illustration, we also plot the result in Fig. [B] 
Again in full qualitative agreement with the analysis in 
SecHH the AMRs in both models are zero for a = 0. They 
also vanish in the limit of a — > oo since no AMR occurs 
if the system is not magnetically polarized. For interme- 
diate ratios of the strengths of the electric and magnetic 
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Magnetization 
of scatterers 


a / 
13 = 


a = 
/3/0 


a = 13 


Along 


[100] 


( aQ-%A \ 


^ _ / cro \ 
V cro - §5 j 


a = 


/ (^0 \ 
V ao 1 


Along 


[010] 


^ _ / '^0 ) 
I a„ - |A i 


5- _ ( '^o - |B \ 
I ao J 


(T = 


1 ao \ 
V ao i 


Along 


[110] 






a = 


/ ao \ 

V (70 / 


Along 


[110] 






(T = 


( (TO \ 

V (70 / 



TABLE I: Conductivity tensor for a 2DEG confined in a [001]-grown III-V semiconductor heterostructure at different 
magnetization of scatterers. Here, ao = e^neT/m, see also Eq. (|A7p . Tie is the electron density, see Eq. (jASp . and 
Ajc? = B / (3'^ = e^uiT /{■kK'^). The conductivity corrections depend essentially on the type of spin-orbit interactions which is 
either Rashba {a) or Dresselhaus (/?) one. The conductivity expressions for arbitrary a and (3 can be found in Appendix 1X1 



parts of the potential, a positive AMR in the Rashba 
model reflects the tangential spin- 1/2 texture while the 
negative AMR in the Dresselhaus model reflects the ra- 
dial texture of the states with large group velocity pro- 
jection to the direction of the current. The singular peak 
at a = 1 originates from the coherent superposition of 
non-magnetic and magnetic scattering amplitudes which 
results in zero scattering probability of one of the two 
states moving along the current direction;^ as we already 
pointed out in Sec IIIl and illustrated in Fig.[T]Jc). 



for Rashba and Dresselhaus models and the magnetic po- 
tential with M oriented along the main in-plane crystal 
axes and along the in-plane diagonals (derived as shown 
below). The component an in the table corresponds to 
the longitudinal response to E along the [100]-axis and 
CT22 along the [010]-axis. To obtain AMR values for elec- 
tric field along an arbitrary angle (j) measured from the 
[100]-axis the conductivity tensors with appropriate mag- 
netization direction of scatterers have to be rotated by 
R-^aR^ where the rotation matrix is given by, 



Rd, — 



[a] 





(b) 




1^2 3 



FIG. 6: AMR for current flowing along [100] crystal axis in 
a pure (a) Rashba and (b) Dresselhaus systems with electro- 
magnetic impurity (oc al + a^.y), varying a, the ratio between 
the electric and magnetic part of the potential. The Fermi en- 
ergy Ef is taken much larger than the spin-orbit interaction, 
so that the Fermi radii of the two bands become almost equal, 
see details in text. 

In Tab. HI we included conductivity components ob- 
tained from the exact solution to the Boltzmann equation 




(17) 



The AMR as defined in Eq. ^ is independent of (p in the 
Rashba model confirming the absence of crystalline AMR 
components in this system. In the Dresselhaus model, 
AMRs of opposite sign are obtained for current along the 
main in-plane axes {4> = 0,7r/2) and along the diagonals 
{(f> = 7r/4, 37r/4), consistent with the crystalline nature of 
the AMR inferred in Sec |TT1 A closer inspection of the 
full angular dependence of the AMR in the Rashba and 
Dresselhaus models allows us to relate our quantitative 
microscopic results to the standard phenomenology of the 
angle-dependent longitudinal resistivity for systems with 
cubic anisotropics,— 



P{'^,(f')/Pav - 1 



Ci cos2{uj ~ (j>) + (18) 
Ci^c COS 2{uj + (/>) + Cc cos 4w , 



where uj and (j> denote the direction angles of M and E to 
the [100] crystal axis, respectively, and pav is the average 
resistivity over all magnetization directions. The coeffi- 
cient Ci of the non-crystalline AMR component, which 
depends only on the relative angle between current and 
magnetization, equals 1/3 for the Rashba model and 
for the Dresselhaus model. The coefficient C/_c of the 
first crystalline component is non-zero (equals —1/3) in 
the Dresselhaus model and zero in the Rashba model, 
consistent with the crystalline nature of the AMR in the 
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[100] 1 + ii' cos(7fc + 7fcO 

[010] l-M'cos(7fc + 7feO 

[110] l-ii'sm(7fe + 7fe/) 

[110] 1 + ii' sm(7fc + 7fcO 

[001] 1 - m' cos(7fc - 7feO 

TABLE II: Magnetization-direction-dependent factors P"''' of 
Eq. (|21|l relevant for magnetic impurities. Functions cos7fc, 
sin7fe are given in the main text. 

Drcsselhaus SOI system and non-crystalline AMR of the 
Rashba system. The coefficient Cc of the higher order 
crystalline term is zero in both models. 

We conclude this Section by presenting the exact so- 
lution to the Boltzmann equation and the corresponding 
AMR values for the combined Rashba-Dresselhaus model 
which implies the dispersion law i?±,k = i kng, 

where Kg = A/a2"T"^M-"2a/3siii2^ is the ^-dependent 
subband spin splitting. We consider a general case of 
arbitrary a and (3 but restrict ourselves to the pure mag- 
netic impurity potential. The electron group velocity 
(l/?i)Vki?k± is now anisotropic and given by 

v±,fc|a- = hkr,/m± (/3cos7fc -f asin7fc)/fi, (19) 
v±,fc|j^ = hky/m± (acos7fc + f3sin^k)/fi (20) 



with kx = kcosO, ky = ksinO, and sin7fc = {acos0 + 
f3 sin 9)/ Kg, cos 7^ = {(3 cos 9 + a sin 9)/ Kg. 

The derivation relies on vanishing angular integrals of 
the generating functions of z/;(i, k; i', k') cx P{i,9;i' ,9') 
(summarized in Tab. |TT] for M along the main in- 
planc crystal axes and the in-plane diagonals) which are 
cos 9 /Kg and sin 9 /Kg. As in the case of the Rashba model 
discussed above, the independence of ^ d9'w{i, 9; i', 9') 
on 9 implies that the non-equilibrium distribution func- 
tion contains only the group velocity, see Eqs. p9l20p . 
and the generating functions of P{i,9;i' ,9') which are 
cos 9 /Kg and sin9/Kg. Note that for arbitrary a and /3 
and for the orientations of M considered in Tabs. |T] and 
|TT]the transition probabilities can then be written as. 



w{i, k; i'k') = —P^^' (i, 9; i'9') , (21) 



where v — ra/Tih^ is the density of states, the k- vector 
independent constant r is given by Eq. (jlip . and the 
angular probabilities P"^^' {i, 9; i'9') are explicitly written 
in Tab. [ill The integral Boltzmann equation ([8]) is then 
solved by the distribution function of a form 



/(z,k)-/°(£;,.k) = T|e|E.v,,k 



dE, 



j,k 



r|e| df{E,M) 



dE, 



j,k 



.Cm 



COS f 
Kg 



Kg 



Ex, 



Kg 



, ^ sm t , _ 
" Kg 



(22) 



Values of the coefficients a%^''y , h%'-''y depend on the magne- 
tization vector direction cm and are given in AppendixIXl 
For jaj = |/9|, and general a, /?, analytical expressions 
for the conductivity tensor of the Rashba-Dresselhaus 
model and short-range magnetic impurity potential with 
M oriented along the main and diagonal in-plane axes 
can be found in Tab. HI and Tab. IIIII in Appendix, re- 
spectively. As pointed out in Sec. [ill the AMR vanishes 
for jaj = \(3\. For a 7^ /3, however, the AMR is non-zero 
and depends both on the relative angle between current 
and magnetization and on the direction of current with 
respect to the crystallographic axes. The AMRs for var- 
ious current directions can again be calculated by rotat- 
ing the conductivity tensor given in Tab. IIIII For current 
along the [lOOj-axis, e.g., and jaj > we obtain 

2(1 _ r2)2 

^^-^ " 2(1 + r2)2 + (3 + r^)n^EF/{ma'^) ' ^^^^ 
where r = fi/a. In the opposite case of jaj < the 



result is the same up to an exchange of a and (3 in Eq. (j23p 
and in the definition of r. 

The smooth transition of the AMR from the pure 
Rashba to pure Drcsselhaus model described by Eq. (|^51) 
is shown in Fig. [71 for Ep = and for intermediate Ep 
corresponding to both majority and minority Rashba- 
Dresselhaus bands occupied. We point out that for a 7^ 
the AMR originates from not only the anisotropic spin- 
texture on the Fermi contours but also, unlike the pure 
Rashba or pure Drcsselhaus models, from anisotropic 
group velocities. In the special case of jaj ~ \(3\, these 
two sources of anisotropy disappear and AMR vanishes 
for any short-range electro-magnetic potential. 

The relative displacement along the diagonal direction 
of the two circular Fermi contours is nevertheless a sig- 
nificant remaining imprint of the SOI in the band struc- 
ture of the jaj = \(3\ model. The AMR can reappear 
if w(«,k;i',k') picks up a dependence on k and k' due 
to other than the spin-texture effect. As pointed out in 
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° p/a a/p ° ° p/a a/p ^ 

(a) (b) 

FIG. 7: AMR for pure magnetic potential impurity as a func- 
tion of the ratio a/f3. Pure Rashba (Dresselhaus) interaction 
corresponds to the left (right) edge, (a) Single band case, 
Ef = (very low electron concentration). Insets show the 
spin textures for several chosen values of q//3. (b) Two band 
case with EpjoP" > fixed. Fermi lines are shown schemat- 
ically, spin textures of the majority band are qualitatively 
similar to the single band case. In the limit Ep oo, the 
AMR vanishes for any value of a/p. 

Sec. ini a long-range (electro-)magnetic impurity poten- 
tial combined with the two displaced Fermi circles would 
yield wave vector dependent w(i, k; i', k') and a non-zero 
AMR even for \a\ = |/3|. 

IV. DISCUSSION 

Calculations in the previous sections show the follow- 
ing trends in the AMR: (i) For the Rashba-Dresselhaus 
model with a short-range magnetic impurity potential, 
the AMR is large (100%) when the minority band is 
depleted and when the SOI is of a pure Rashba type 
{(3 = 0) or pure Dresselhaus type (a = 0). (ii) The 
AMR vanishes when |a| = |/3| or for an arbitrary a and 
/3 when the majority and minority bands become nearly 
degenerate, (iii) For impurities containing a combined 
electro-magnetic potential, the AMR has the same sign 
as for the pure magnetic impurity potential, is maximized 
when the two components have equal strength, and re- 
mains large (200%) even in the limit of nearly degenerate 
Rashba-Dresselhaus bands, (iv) We have also noted (in 
agreement with Ref. [29| ) that in the higher-spin Kohn- 
Luttinger model, the AMR is expected to have oppo- 
site signs for pure magnetic potential and for electro- 
magnetic potential with comparable strength of the two 
components. Wc will now discuss implications of obser- 
vations (i-iii) and inspect the applicability of our linear- 
response quasiclassical theory for 2D systems with real- 
istic material parameters. 

Two-dimensional electron systems with Rashba and 
Dresselhaus SOI have been studied in n-type InAs and 
GaAs quantum wells^Si'^i^i^ with mobilities iJ. up to 
3 X 10^ cm^/Vs and 3.5 x 10^ cm^/Vs, and magnitudes of 
the SOI of the order of - 10"^^ cVm and - 10"^^ ^Vm, 
respectively. The ratio |a//3| is ranging between approx- 
imately 1.5 to 8 for these two-dimensional systems with 



electron densities of the order of ~ 10^^ — 10^^ cm~^. 

The semiclassical Boltzmann theory is applicable when 
the following two conditions are satisfied. First, the par- 
ticle's de Broglie wavelength must be smaller than the 
mean free path. At low temperatures (as compared to 
the Fermi temperature) the condition implies that 

> ^ , (24) 

where r = m/i/e. For the above InAs and GaAs 
two-dimensional systems^ m/{hT) is of the order of 
lO^*' cm~^ and lO^cm"^, respectively, so the inequality 
p4)) can be safely met. 

The second condition requires that the smearing of the 
spin-split bands due to disorder is smaller than the spin 
splitting energy i?+,k — -E^-,k- Since the AMR we study 
is due to the SOI in the band structure (rather than 
in the scatterers) it remains non-zero only in the strong 
SOI/weak disorder regime. As a consequence, the con- 
centration must also fulfill the following inequality 

rie > H'^/SttkIt'^. (25) 

Assuming a pure Rashba system (i. e. Kg = a), the 
right-hand side in (|25p is of the order of 10^ cm^^ for 
both InAs and GaAs, respectively, so the condition is 
again satisfied for typical electron densities. Introduc- 
ing magnetic impurities will certainly decrease the mo- 
bility of the two-dimensional systems, nevertheless, con- 
ditions and (|^5|) might remain satisfied for feasible 
electron densities. We also note that the inequality (j^S]) 
can be reformulated in terms of the mean free path Ir and 
spin precession length, which can be roughly estimated 
as Xs ~ f?j{rna). Namely, must be larger than so 
that an electron randomizes its spin orientation due to 
the spin-orbit precession between two subsequent scatter- 
ing events. This restriction corresponds to the approx- 
imation which neglects the off-diagonal elements of the 
non-equilibrium distribution function in the spin space. 

Having established parameter range of the validity of 
the Boltzmann approach wc can now return to points 
(i-iii) from the beginning of this section and comment 
on the expected AMRs for realistic material parameters. 
Since for short range impurities the AMR is weak when 
|a| « 1/3 1 let us assume pure Rashba model only. By a 
direct inspection of the results in Tabs. HITTIl we find that 
the ratio between the isotropic and anisotropic part, gq 
and CTi, of the conductivity tensor depends on the SOI 
strength and electron density and can be estimated as 




For usual electron densities ~ 10^"'^ cm^^, this ratio 
will be of the order of 0.01 for a pure magnetic im- 
purity potential, implying weak AMR of the order of 
1%. By depleting the minority band, the ratio ai/cro 
can be enhanced and the AMR can reach up to 100% 
(recall Fig. [5]). However, corresponding densities of 
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Ue « 10^ cm ^ are relatively low compared to densities we get 



of typical experimental two-dimensional electron systems 
and also we then move towards the edge of the validity 
of the Boltzmann theory. 

The AMR will be further reduced by the presence of 
another impurities than the (electro-)magnetic ones. In 
terms of resistivities, this follows from the Matthiessen's 
rule stating that the total resistivity is a sum of resistivi- 
ties due to the particular scattering mechanisms 4^ Since 
scattering from pure non-magnetic impurities yields zero 
contribution to the type of AMR discussed in this paper 
the overall relative magnetic anisotropy of the resistivity 
is suppressed by their presence. 

On the other hand, for impurities containing a com- 
bined electro-magnetic potentials which add up coher- 
ently during the scattering, the AMR is expected to be 
largely enhanced even in the high density regime. The 
strongest AMR is predicted for similar strength of the 
magnetic and electric parts of the scattering potential. 
This applies, e.g., to Mn in GaAs which acts both as a 
charged dopant and a localized magnetic impurity, and 
the AMRs in GaAs:Mn can reach ~ 10%4i In the present 
paper, we however wish to limit our investigation of mod- 
els beyond the Rashba-Dresselhaus one to the qualitative 
discussion of Sec. Ill Dl complemented by the exact Boltz- 
mann equation AMR given in Appendix |B] Wc refer the 
reader to Refs. [?9ii37..38 for a more quantitative discus- 
sion of AMR in (Ga,Mn)As and finally remark that the 
realization of large AMRs in Rashba-Dresselhaus systems 
with electro-magnetic impurities will require doping with 
magnetic donors. 
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APPENDIX A: BOLTZMANN EQUATION FOR 
RASHBA-DRESSELHAUS HAMILTONIAN 



In order to determine the non-equilibrium distribution 
function, we insert the ansatz (|22p into the Boltzmann 
equation ([8]) and obtain a set of four linear equations for 



parameters 



h%^'y, one for each direction of electric 



field (a;, y) and each magnetization direction ejv/. For the 
scatterers magnetized along i(y)-axis and E = (i?a;,0). 



, 1 

ax = T ^ ± ^ 



"|a2_/32| 



(Al) 



2a/3 |a2_^2| 
(a2 + /32)2 /|a2-/32 | 



2 



4a/3 V a^+P^ 

2 i «2 / ^2^p2 



- 1 



(A2) 



a:r— TT- 1 



2aP \' I a2 - /?2 I 
The choice E ^ (0, Ey) leads to 

1 r p2_^2 

±2 



bx 



a2 -I- /32 



-|- a 



\a^-l3^\ 
a^ + l3^ , 1 



y < 2 



± - 
2 



1 - 



/?2-a2 
^ 2a(3 



I a2 -/?2 |_ 

(A3) 



I a2 _ ^2 I 



(A4) 



^ 2al3 \a^-l3^\. 

Here, we skip the cm superscript for brevity and relate 
the upper and lower signs to the magnetization M along 
X and y axes respectively. 

For scatterers magnetized along the [110] axis we have 

a2+/32- |a2-/32| bx 



a(3 



Aa(3 ' 2 

^ ^ (a2+/32)2 / |a2_^2| ^ 



4q2/32 ^0,2+^3 



(A5) 



bx 



a2+/32- I a2_^2 I 



AaP 2 2 

while equations for ay and hy can be obtained from 
Eqs. (|A5IA6[) by the substitution ax —> by, bx —> ay. 



1. Impurity magnetization along the [100]-axis. 

Here, we assume that impurities are magnetized along 
the X-axis, i. e. the scattering potential is proportional 
to ax- 

li a > (3 then solution of Eqs. (|A1IA2|) and 
Eqs. (|A3IA4[) with upper sign reads 

ax = —, bx = ^ 2 , a2 ~ 2"-"' 



3a2 + /32 



3a2 + (3' 



and 

„ n 8a3/3 , 24^4 
In the opposite case f3 > a the coefficients are 



- 4/32^2 ^ 3^4 ^ 2a(a2 - f3^) 



a2 + 3/32 



2Q(a2 - /?2)/3 /34_Q,4 



a2 + 



^ a2 + 3/32 
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Magnetization 
direction 

[100] 



[010] 



f3>a 



[110] 



[110] 





2e2mc«(c«2-^2j^^ . 




2e^m.0^(a^-f3'^)T 








2&^ma{a^ —i3^)[3r 


fi47r(3Q2+/32) 


fi<l7r(3Q2+/32) 
2e2ma2(a2_/32)^ 


fil7r(3Q2+/32) 


ft-l7r(3Q2+^2) 


e2m(a-/3)(a + /3)2T 


e2m(a-/3)(a + /3)2r 


e2m(a-/3)(a + /3)2T 


fi<l7r(3a + /3) 
e2m(a-/3)(a + /3)2r 


ft47r(3a+/3) 


fi'l7r(3a + /3) 




e2m(a-/3)(a + 0)2T 


fi*7r(3a+/3) 
e2m(ct-/3)(a + /3)2r 


ft"l7r(3c« + /3) 
e2m(a-/3)(a + /3)2r 


ft'»7r(3c« + /3) 


fi<l,r(3a + /3) 



2e2ma2(n2_^2-j^ 
2e2ma(a2-ff2)^^ 



' 2e^m{a^-IS-')fi-'T 
fi4^(jj2_^3^2') 
2e2mct(ct2-;32j^^ 

fi47r(Q2+3/32) 

e2m(a-/3)(a + /3)2T 
fi47r(a + 3/3) 

e2m(Q;-/3)(ct + /3)2T 
fi4-ir(a + 3/3) 



m(a- f3)(a + T 
fi47r(c« + 3/3) 
e2m(a-/3)(ct + )3)2T 
ft4^((i + 3/3) 



2e2ma(a2-/32)^^ ^ 

ft4^(c«2+3/32) 
2e2m(a2-/32)^2^ 

ft'l^(a2+3^2') ^ 

2e2ma(a2-/32-)^^ > 

ft4^(c«2+3/32) 
2e2>na2(a2-;32.j^ 

ft4 7r(a2+3,32) 

e'^m(a-l3)(oi + l3)^T 

e2m(a-/3)(a + /3)2T 
fi47r(a + 3/3) 



e2m(a-/3)(a + /3)2 
fi47r(a + 3/3) 
e^m(a-f3)(a + f3)^T 



TABLE III: Anisotropic part di of the total conductivity tensor a — Icro + cfi for arbitrary a and (3. 



2. Impurity magnetization along the [010]-axis. 

The scattering potential is proportional to ay in this 
case, li a > (3 then the solution of Eqs. (jAllA2p and 
Eqs. (|A3IA4[) with lower signs is given by 



24a^ 



while in the opposite case {P > a), the coefficients read 



a,- = 



a„ = 



h3/32 
a2 + 3/32 



, by 



a2 + 3/32 ' 
- 4/32^2 + 3/?^ 
a2 + 3/32 ■ 



3. Impurity magnetization along the [110]-cixis. 

Here, the scattering potential is proportional to ^{o'x- 
Gy). li a > 13 then the coefficients a,j; y and h^^y read 

- 3/3a2 ^ ^2^ ^ ^3 2a2 (a - /3) 

= , Ot = 

3a + /3 3a + /3 

_2a2(a-/3) _ - 3/3a2 + /32a + /J3 

3« + /3 ' 3^T^ ■ 

In the opposite case /3 > a we have 

+ /3a2 - 3/?2a + , 2/32 _ q,) 

0:7: = :rT. , 0,. — 



a + 3/3 



a + 3/3 ' 



_ 2/32(/3-a) __a3_+/3a2-3/32a + /33 



a + 3/3 



a + 3/3 



The case when the impurities are magnetized along the 
[110]-axis can be treated in the same way. 



To write down the conductivity for arbitrary a and /3 
it is convenient to define 



CTo 



(A7) 



where the electron concentration at Ep > can be ex- 
actly expressed as 



/32 



(A8) 



Thus defined (Tq becomes identical with the Drude for- 
mula when a = /3 = 0. In fact, ctq times unity 2x2 ma- 
trix describes the conductivity of a 2DEG due to the non- 
magnetic short-range scatterers, see Ref. [l3- In the pres- 
ence of magnetized scatterers the conductivity acquires 
an additional term di which is summarized in Tab. IIIII 
To obtain the total conductivity tensor one has to sum 
up both these terms, i.e. a = Ictq + 0^1 • Conductivity 
tensors under special conditions in Tab. [T] can be recov- 
ered by a proper choice of a, (3 in Tab. IIIII Table IIIII 
thus summarizes the main computational results of this 
paper. They describe an additional term in the electrical 
conductivity of a 2DEG confined in a [001]-grown III- 
V semiconductor heterostructure due to the magnetized 
elastic scatterers. 



APPENDIX B: BOLTZMANN EQUATION FOR 
KOHN-LUTTINGER HAMILTONIAN 

Non-trivial exact analytical solutions to the Boltzmann 
equation ([8]) exist also for some models in d = 3 dimen- 
sions. The one described in Section fll DI constitutes one 
such example and we outline here the main steps needed 
to calculate the AMR in this model and thus confirm the 
appropriateness of the sketch on Fig. IH^b) . 

Physically, the model concerns carriers of the two 
heavy-hole Eg bands (HH bands) scattered off mag- 
netic impurities. The band structure can formally be 
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viewed as the (71 — 272)7(71 + 272) limit (negligible 
light-hole density of states) of the Hamiltonian ^ with 
jx,y,z = Ssx.y^z while scatterers uniformly polarized along 
z-direction are modeled by y oc Sz in terms of Eq. 
Explicit expressions for the spin matrices Sx,y.z can be 
found e.g. in the Appendix of Ref. |4^. Without going 
into details, we remark that this model could be used 
to describe the AMR in (metallic) p-type III-V or II- VI 
semiconductors with dilute Mn impurities if their charge 
is either zero or strongly screened;— i^i^ Mn atom d- 
states hybridize with the host valence band and create^ 
a fc-independenti^ impurity potential V oc em • s. Rele- 
vant values of the proportionality constant and host ma- 
terial band structure parameters 71 , 72 can be found in 
Ref. [13. We also stress that we will be treating a model 
where the densities of states of the two involved (HH) 
bands are equal {h 0) and this is of course (again) 
only an approximation to realistic systems. 



Non-equilibrium distributions due to applied electric 
field turn out to be the same for both HH bands in such 
a model, and we are required to solve three decoupled 
integral equations 



^47r/3TYi"(r2) = w{n)c{n)- J dn' w{n,n')c{n') 
-^/Sn/STY^^in) = w{n)p{n) - J dn' w{n,n')p{n') 
y^8n/3TYf^{n) = w{n)q{n) - J dfi' w{n,n')q{n') 

(Bl) 

where r is defined by the 3D analogy of Eq. (fTTjl . de- 
notes a compound variable tp, •& parameterizing the unit 
sphere, so that J dfl — sin ?9 d-d J^^ dip, 
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and Y,™ = = Yri^') denote spherical harmonics normalized to / dl7r/"'*(f7)r;™'(f7) = dwSr, 

according to the Condon-Shortley convention. 

Non-equilibrium distribution under the effect of E = {E^, Ey, E^) is then 



(B3) 



for the both bands (which arc in the h ^ approxima- 
tion identical) , and conductivities implied in the spirit of 
Eq. HH) are 

dn 1 {p{n) + q{n)) (rf ^ (n) - y^ {n)) 

'yy - , df^^(p(f^)-9(f^))(Yi^(f^)+rfi(f7)) 
dVL^c(Sl)Y^{^) . 



in units of ^Ji|2n ■ e^rie/m. Using analytical solutions of 
Eqs. (jBip corresponding to scattering amplitudes (|B2[) . 
we find 



AMR =2^^±^^^^i:2^« 0.45, 
2 + V2 arctan ^/2 



(B4) 



according to our definition of AMR ([T]), that is resis- 
tance parallel to the magnetization is higher. It is thus 
confirmed that sketches for pure magnetic scattering 



in Fig. Ub) appropriately describe conductivity calcu- 
lated by exactly solving the Boltzmann equation. We 
also obtained the same result (jB4|) within the Keldysh 
formalism2ii2^ where conductivity Uxx turns out to be 
proportional to {vx5vx5)t with Vx ~ dHuL/dkx, S = 
qr _ qA^ qA ^ ^QRy ^j^j (Gfl)-i =Ef~ Hkl - 

with self-energy (x s^. Brackets (. . .)t in this Kubo- 
formula-type result^ mean trace in the space of 4 x 4 ma- 
trices and integration over the fc-spacc. Conductivity cr^z 
is analogous {vx is replaced by Vz, S'^ stays unchanged). 



We finally remark that Eqs. (jBTj) are completely analo- 
gous to the two equations (8,7) of the 2D case in Ref. [3^. 
Solution of those equations was constructed in the form of 
a Fourier series or a modified Fourier scries as explained 
in the note added in proof of that reference. In our cur- 
rent 3D problem defined by Eqs. (jBllB2p . if expanded 
in terms of modified spherical harmonics Yi"^{n)/w{n), 
the solutions c{il), p{fl), q{fl) are found to contain only 
few terms [the harmonics present in the left-hand-sides of 
Eqs. (jBip : note the analogy to the discussion of Sec. lIII Al 
and Eq. (|22p when spherical harmonics replace sines and 
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cosines]. During the calculations we have to be how- simple scalar products of spherical harmonics where or- 
ever cautious as the integrals of Eqs. (|B1[) do not contain thogonality relations apply. 
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